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Methodology for the study of modified jet-like topologies in heavy ion collisions 

via three particle correlation functions 
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Methodology is presented for analysis of three-particle correlation functions obtained in heavy ion 
collisions at ultra-relativistic energies. We show that harmonic correlations can be removed and jet 
driven correlations reliably extracted. Results from detailed Monte Carlo simulations are used to 
demonstrate the efficacy of this technique for the study of modifications to away-side jet topologies. 
Such modifications are an essential probe of the properties of the quark gluon plasma produced in 
heavy ion collisions. 

PACS numbers: PACS 25.75.Ld 



I. INTRODUCTION 

Studies of ultra-relativistic heavy ion collisions have 
provided a wealth of evidence for the creation of a new 
state of matter in a collision zone of very high energy 
density [l], H, H H| , historically termed the quark gluon 
plasma QGP. As the properties of this matter are being 
explored and characterized, it may be that new termi- 
nology will be chosen to be more suggestive of its actual 
properties, as they are revealed. Current research is ac- 
tively directed toward measurements of these properties. 

One of the major tools for such measurements has 
been the study of correlations between the observed par- 
ticles which emanate from the collision medium following 
its expansion and ultimate hadronization. For example, 
much detailed information has been obtained for the so 
called "harmonic flow" correlations between individually 
selected particles and the reaction plane (e.g. the second 
harmonic coefficient (v^)) as a function of particle iden- 
tity PID, collision centrality, transverse momentum px, 
and rapidity 77 g & & S M M 110- From 
the set of systematic data, initial estimates have been 
made of the speed of sound in the QGP, the ratio of its 
viscosity to entropy density, its bulk viscosity and other 
properties 0, M M El M M M M 11 • 

Formation of the collision medium is sometimes accom- 
panied by hard parton-parton scatterings. These scat- 
tered partons can interact strongly with the medium and 
lose energy as they propagate thr oug h it, before frag- 
menting into jets of hadrons [23|,[24],[25[|. Such energy loss 
can lead to a strong modification of both the yield and 
the topological patterns of jets [H |27|. This gives an- 
other detailed experimental probe of the medium which 
is currently being explored and developed via measure- 
ments of di-jet correlations [H, US HI [H H [H| . 

Techniques for the study of jet-induced two-particle 
azimuthal angle correlations are well advanced, and re- 
sults are being accumulated for systematic evaluation. 
In p+p and d+Au collisions, these jets of hadrons are 
found essentially back-to-back in azimuth i.e. A<fi ~ 
180° [32} . By contrast, the observed di-jet topologies 
in Au+Au collisions are found to be significantly mod- 



ified [2a I23, |28|, 12a |3fJ, l3l|, presumably due to parton 
medium interactions that occur on passage of the parton 
through the reaction medium. One very intriguing result 
of these di-jet studies is the observation of broadening of 
the away-side jet, and even the displacement of its most 
probable angle away from = 180° [H[27]]. Several 
mechanistic scenarios have been proposed for this obser- 
vation; they include Cerenkov gluon radiation [33 j. coni- 
cal flow [Mil [Mil [Mil and deflected or "bent" jets 
[40L l4lj . To date, the characteristic px dependent pat- 
terns predicted for Cerenkov gluon radiation [33[ have 
not been observed. The confirmation of a conical flow 
signal would not only give a direct probe of the equation 
of state (EOS) of hot QCD matter, but also an impor- 
tant constraint for an upper limit for the viscosity of the 
medium [H [Hil i|. 

Two particle correlation measurements do not provide 
an unambiguous distinction between conical flow and de- 
flected jets. However, the topological information af- 
forded by the correlations between three or more particles 
can. Here, we lay out a method of analysis for three parti- 
cle correlations which demonstrates a topological distinc- 
tion between conical flow and deflected (or bent) jets. We 
follow and build on methodology formerly presented for 
the study of two-particle correlations; namely, the use 
of normalized correlation functions and extensive test- 
ing via Monte Carlo reaction simulations [42]. To focus 
on the di-jet-like characteristics, we exploit a novel co- 
ordinate system that is most intuitive and natural for 
visualizing di-jet topologies. 



II. TWO-PARTICLE CORRELATION 
FUNCTIONS 

In reference (42[ we presented a method for construct- 
ing and analyzing two particle correlation functions based 
on the relative laboratory azimuthal angle A(f>, for par- 
ticle pairs. In brief, jet correlations were emphasized by 
selecting events with at least one high transverse momen- 
tum pt (trigger) particle. Each trigger particle was then 
paired with associated particles of a lower px to obtain 
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the pair correlation function C2(A</>); 



C 2 (A4>) 



(1) 



where iVg(A0) and Nm{A4>) are normalized same-events 
and mixed-events distributions and A(f> = \<p± — fa] is the 
difference between the azimuthal angles of the particle 
pair. The same-events distribution was constructed from 
particle pairs obtained from the same event; the mixed- 
events distribution was constructed by selecting each par- 
ticle in a given pair from a different events having similar 
centrality and collision vertex positions. 



A. Extraction of Jet Shapes 

Using a two component model ansatz, we showed [42[ 
that the pair correlation from a combination of flow and 
jet sources is given by; 



C 2 (A0) = & o [C//(A0) + Cj{A<j>)], 



(2) 



where the flow contribution Ch(A4>) can be estimated as 

C H (A<t>) = [1 + 2vlcos2{A(j)) + 2v\cosA{A(j))}, (3) 

and Cj(A</>) is the jet function that needs to be eval- 
uated. It is noteworthy that no explicit or implicit as- 
sumption is made for the functional form of Cj(A(f>). By 
rearrangement of Eq. [2] one obtains 



Cj(A^) 



C 2 (A0)-6 o C g (A0) 
6o 



(4) 



Thus, knowledge of bo is required to evaluate the jet-like 
function Cj(A<t>). It is clear from Eq. 2] that bo is in- 
fluenced by the jet function which is being sought after, 
and an approximation is required to aid the evaluation of 
bo. After extensive detailed simulation studies, we con- 
cluded that a reasonable assumption for the extraction 
of reliable jet-like functions is the zero yield at minimum 
(ZYAM) condition - that the di-jet function has a zero 
yield at minimum (42j; 



& Cff(A</> r , 



C2(A</> min ), 



(5) 



which can be solved to obtain b . This procedure allows 
reliable extraction of jet-like topologies and a lower limit 
for jet-like per trigger yields. 



B. Simulations 

Reaction simulations were used to extensively test this 
analysis approach. The requisite simulations were car- 
ried out on an event-by-event basis with the following 
essential steps: 

• First, the reaction plane orientation was chosen. 



• Particles were then emitted with pt and multiplic- 
ity according to the observed distributions. 

• For flowing particles, the azimuthal angle for each 
particle <fii was chosen to give a harmonic distribu- 
tion with respect to the angle of the reaction plane 



N(<P ~ tp R ) oc [1 + 2v 2 cos2{cj) - tp R ) 

+ 2V4C0S4(4> — IpR)} 



(0) 



where «2,4 are Fourier coefficients which character- 
ize the strength of the flow. 

• For jets, the orientation of the lead- or near-side 
axis was chosen with a random azimuth. For "nor- 
mal" jets, the away-side axis was oriented opposite 
to the lead-axis on average. 

• The lead-jet particles were emitted clustered about 
the lead-axis. The away-side jet particles were 
emitted clustered about the away-side axis. 

Simulations were performed both for an ideal detector 
and for the PHENIX detector. A detailed description 
of the latter can be found in Ref. 43]. As outlined 
in Ref. (42| . the simulations indicated that, even for 
cases in which strongly distorted away-side jets were in- 
troduced, the decomposition method retrieved the shape 
of the input jet function in detail, confirming that the 
decomposition procedure is robust. 

We now turn to the discussion of three particle corre- 
lation functions. 



III. THREE-PARTICLE CORRELATION 
FUNCTIONS 

Analogous to C2(A</>), three particle correlation func- 
tions can be constructed from particle triplets comprised 
of one high pt trigger particle (i) and two associated par- 
ticles of lower pr (j and k). To focus on three-particle jet- 
like correlations, we transform the lab angles {Oiab, 4>iab) 
of each particle to (9* , 4>*) in a new frame whose z axis is 
the direction of the high px trigger hadron. This coordi- 
nate frame is illustrated in Fig. [Ha). In this frame, the 
three particle correlation function is given by the ratio of 
two distributions: 



C 3 (0*,A<£*) 



A s (r,A^) 

N M (6*,A<p*) 



(7) 



where N s (6*, A0*) and N M (0*,A^*) are normalized 2D- 
distributions for the same- and mixed-events respectively. 
Here, 9* is the polar angle of one of the two associated 
hadrons and A0* = 1 0* — <\>* k | is the difference between 
their azimuthal angles. The same-events distribution was 
obtained via event-by-event selection of particle triplets 
from the same event. For mixed-events, particle triplets 
were obtained by selecting each member from a different 
event. 
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Away-side Jet 



FIG. 1: (a) Schematic illustration of the coordinate system used for the construction of three-particle correlation functions 
comprised of high px trigger particle and two low px associated particles. The lab angles (8i a b, (j>iab) of each particle is 
transformed to (0*,<j>*) in a new frame whose z axis is the direction of the high px trigger, (b) Polar plot of a simulated 
three-particle correlation function for a normal di-jet. The polar angle 6* , of one of the associated low px particles in the new 
frame, is plotted along the radial axis; the difference between the azimuthal angles A(f>* , of the two associated hadrons in the 
new frame (see text) is plotted along the azimuthal axis. 



The correlation function so obtained is best viewed in 
the polar representation shown in Fig. [TJb) , where a sim- 
ulated three particle correlation surface is shown. Both 
flow and di-jet correlations were incorporate in the sim- 
ulation. In Fig. Wlb), 9* and A<p* indicate the radial 
and azimuthal axes respectively. In this polar represen- 
tation, the near-side jet is indicated by a peak at the 
center of the plot (8* = 0°) and the characteristic ridge 
at 6* = 180° signals a normal or unmodified away-side 
jet (i.e. a back-to-back jet). 

The primary objective of our study is to use such cor- 
relation surfaces to distinguish between different mecha- 
nistic scenarios for away-side jet modification which can 
not be discerned via two-particle correlation functions. 
To demonstrate this ability we used our simulation code 
to model di-jets with (i) an away-side bent jet and (ii) an 
away-side "cone" jet (see illustrations in Figs. |2](a) and 
(c) respectively). For the first, the away-side jet axis is 
bent to an angle of ~ 120° with respect to the lead-jet 
axis and the away-side jet particles are emitted clustered 
about this axis. The axis for the away-side cone jet was 
chosen opposite (on average) to the lead-axis, and its as- 
sociated jet-like particles were emitted so as to mimic a 
Mach cone with Mach angle 9 m ~ 60°. It is important 
to note here that these simulations were performed for 
the PHENIX detector acceptance. Equally important is 
the fact that model parameters for the simulations were 
tuned (with insight from experimental data) to give the 
same shape for the simulated jet-pair correlation func- 
tions for bent- and cone jets, as shown in Fig. [DJb). For 
this simulation set, an isotropic underlying event was em- 
ployed. 

Figure [Ha) gives an illustration of a di-jet with a bent 
away-side jet. The corresponding simulated three par- 



ticle correlation function is shown in Fig. [DJd). This 
correlation surface exhibits a sizable peak at 9* = 0° 
corresponding to the lead- or same-side jet, and a ridge 
at 9* — 120° corresponding to the away-side jet, shifted, 
on average, by ~ 60° from 9* — 180°. The projection 
of A(f>* for 9* ~ 120° (i.e around the ridge) is shown in 
Fig. [2je); it shows a relatively large peak near A<fi* = 0°, 
which results from particle triplets comprised of a high 
px trigger from the near-side jet and two low px associ- 
ated particles from the away-side jet. Two small peaks 
can also be observed at A0* - 90° and 270° in Fig.^d); 
they result from particle triplets in which the high px 
trigger and one low px particle is from the same-side jet, 
but the second low px particle is from the away-side jet. 
The fall-off at large Acf>* angles reflects the influence of 
the PHENIX detector acceptance as discussed below. 

Figure [2]Jf) shows the three particle correlation surface 
for a di-jet with an away-side Mach cone (see illustration 
in Fig. [He)). Analogous to Fig. [DJd), a near-side peak 
at 9* = 0°, and an away-side ridge near 9* = 120° is ap- 
parent. However, in contrast to the case for the bent jet, 
the peak at or near Acf>* = 0° is less pronounced. This 
distinctive signal is made more transparent in Fig. (2Jc) 
where the azimuthal projections (Atfi* along the ridge) 
for the away-side bent jet and Mach cone are compared. 
These correlations result from particle triplets comprised 
of a high px trigger from the same-side jet and two low px 
particles from the away-side Mach cone. As in the case of 
the bent jet, the two small peaks at Acf>* ~ 90° and 270° 
in Fig. [3ff ) , reflects "anomalous" particle triplets in 
which the high px trigger and one low px particle is from 
the same-side jet and the second low px particle is from 
the away-side Mach cone. 

Intuitively, correlated triplets which result from a near- 
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Deflected jet v '< rad > 

Mach Cone 




FIG. 2: (a) Schematic illustration of a bent jet topology; (b) comparison of the two-particle jet functions for bent jet and 
conical flow; (c) schematic illustration of conical flow topology; (d) simulated three particle correlation surface for bent jet; 
(f) simulated three particle correlation surface for conical flow. Panel (e) shows an azimuthal projection of both surfaces for 
6* ~ 120° i.e along the ridge (see text). 



side high px trigger and two associated particles from an 
away-side Mach cone, should lead to an away-side az- 
imuthal ridge at 8* — 180° — 9m with no preference for 
Acf>* . Therefore, the correlation surface in Fig. [UJf), as 
well as its A<p* projection in Fig. [^e) gives an indication 
of the relative influence of the PHENIX detector accep- 
tance (especially in the vicinity of (j>* angles close to 180°) 
and the anomalous triplet correlations discussed above. 



A. Suppression of harmonic correlations 

The correlation functions shown in Fig. [5] were ob- 
tained from simulated events with hits recorded in the 
PHENIX acceptance from three di-jet-like particles, one 
trigger (high pr) and two associated (lower pr) parti- 
cles; the underlying event particles were made isotropic. 
It is well known however, that there are strong harmonic 
correlations between the various particles and the reac- 
tion plane that follow Eq. [6] Therefore, the underlying 
event is not actually isotropic. Detailed systematic mea- 
surements have been made of the Fourier coefficients ^2,4 
with a significantly reduced influence from jet-like corre- 
lations [1J, |44j; thus, it is straightforward to incorporate 
these harmonic correlations into the event simulations. 

A three particle correlation surface which results from 
the combined influence of jet-like and flow correlations 
is shown in Fig. [31a). It shows that the combined cor- 
relations make it more difficult to visualize the detailed 



effects of the di-jet-like emissions in the raw correlation 
function. This is akin to the actual experimental situa- 
tion. Therefore, we wish to remove the harmonic corre- 
lation effects from the three-particle correlation surfaces. 
To this end we follow a procedure similar to that de- 
scribed in reference [42j using the ZYAM assumption i.e. 
we assume that the jet correlation function has zero yield 
at its minimum. 

The three-particle correlation function which results 
from a combination of flow and jet-like sources can be 
given as 

C 3 (6*,Acj>*) = a [C H [6* , A0*) + Cj{6* , A0*)] . (8) 

The harmonic flow correlation Ch(0* , A(f>*) is obtained 
from a flow only simulation where particles are emitted 
according to the measured pattern with respect to the 
reaction plane. The jet correlation is then given by 

C J (0\Ar) = Cs(e ^ Ar) - aoCH{e ^ ) . (9) 

a 

Applying the ZYAM condition one obtains 

a o C H (0* min , A4>* min ) = C-s{6* min , A(f>* min ) (10) 

which can be solved to obtain a . 

Tests for flow removal for the case of three-particle cor- 
relation functions are shown in Fig. [3] Fig. G^a) shows 
the correlation surface obtained from a simulation that 
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FIG. 3: (a) Simulated three particle correlation surfaces; the simulation included both flow and jet-like correlations, (b) Same 
as (a) but after flow subtraction, (c) Input jet-like correlations only. Note the similarity between (b) and (c) . 



includes both di-jet and harmonic correlations. Fig. [3jb) 
shows the jet function that results after removal of the 
harmonic correlation by the ZYAM-driven subtraction. 
This is to be compared to the input jet-like correlation 
function generated with no harmonic correlations, shown 
in Fig. [3l^c). The similarity between Figs. [U[b) and (c) 
clearly show that the ZYAM procedure is able to suc- 
cessfully recover the essential characteristics of the input 
di-jet-like function. 



B. Removal of (2+1) correlations 

The three particle correlation function Cj (6* , A0* ) , 
obtained after removal of flow effects, still contain con- 
tributions from false triplets or (2+l)-correlations. They 
are of two types: (i) False "hard-soft" (hs) triplets in 
which the high px trigger and one associated low px par- 
ticle come from the di-jet, but the second low px asso- 
ciated particle is from the underlying event, (ii) False 
"soft-soft" (ss) triplets in which the two low px associ- 
ated particles belong to the same di-jet but the high px 
trigger is from the underlying event. An estimate of such 
correlations can be made and removed as follows. First, 
a (2+l)-correlation function 

C (21) (9 ,A0 ) - , (11) 



is constructed from event pairs. Here, the distribution 
for fake-triplets Ns(2+i)(9* , A<f>*) is obtained by taking 
two particles from one event (hs or ss) and a third from 
another event. The hs and ss triplets are sampled in 
the ratio of observed ss and hs two-particle correla- 
tion strengths. As before, the mixed-events distribution 
Nm(0* , A(j>*) is obtained by taking each member of a 
particle triplet from a different event. This correlation 
function estimates both the ss and hs components of 
the (2+l)-contribution. Second, the flow contribution 

was subtracted via the procedure 
outlined earlier, to obtain the (2+l)-jet-like contribution 
Cj( 2 +i)(0*, A(j>*). The fully corrected triplet correlation 
function C$j(9* , A(j>*) was obtained via flow and (2+1)- 
subtraction; 



C 3J (9*,A<t>*) = C 3 ( 



,A<f > *)-H(9* ) A<l>*) 
-C J(2+1) (9*,A<f>*). (12) 



N M (9*,Acj>* 



We have made extensive tests of this analysis pro- 
cedure via detailed simulations and have confirmed its 
utility in suppressing both flow and (2+l)-correlations. 
Fig. 0Ja) shows a representative comparison of the A(j>* 
projections [along the ridge] for an input (filled cir- 
cles) and a recovered (filled squares) three-particle cor- 
relation function. The recovered correlation function is 
obtained after removing the effects of flow and (2+1)- 
contributions. The comparison clearly speaks to the ef- 
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FIG. 4: (a) Comparison of the <f>* projections for 9* ~ 120° for simulated input and output correlation functions. The 
output correlation function is obtained after harmonic and (2+1) subtraction (see text); (b) correlation surface obtained from 
a simulation in which only (2+l)-correlations are introduced i.e no genuine three particle correlations are introduced. A small 
offset is added for clarity. 



ficacy of the technique. A further confirmation of the 
efficacy of our procedure can be illustrated via the corre- 
lation surfaces obtained from simulations which include 
only (2+l)-correlations. One such example is given in 
Fig. HJb); it shows the expected "flat" and featureless 
surface that is to be expected from the successful removal 
of (2+l)-contributions. 

IV. SUMMARY 

Detailed experimental probes of the hot and dense 
medium created in ultra-relativistic heavy ion collisions, 
are currently being explored via measurements of multi- 
hadron correlations at RHIC. These correlations contain 
contributions from both harmonic flow and jet-like emis- 
sions which must be disentangled, as well as topological 
features important to the study of the mechanism for 
modification of the away-side jet. A method is presented 



for the analysis of such data to retrieve strongly modified 
jet-like topologies via three particle correlation functions. 
In particular, we have shown that a distinction between 
conical flow and a deflected (or bent) jet can be obtained 
in our analysis framework. Intuitively, one knows that a 
normal or back-to-back jet would generate a peak for the 
near-side jet centered at 9* = 0° and an azimuthal ridge 
at 9* = 180°. Further, it is expected that a bent jet would 
lead to an away-side azimuthal ridge shifted to a 9* angle 
less than 180°, but would cluster the associated away-side 
jet-like particles near to A<fi* = 0°. By contrast, a Mach 
cone should lead to an away-side azimuthal ridge also 
shifted to a 9* angle less than 180°, but without pref- 
erence for any A0* value. For the most part, these are 
the patterns shown by the simulations, albeit with some 
distortions due to the limited 77 acceptance of the detec- 
tor and the detection of anomalous triplets. However, 
these distortions do not blur the topological distinction 
between an away-side bent jet and conical flow. 
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